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ONE-DIMENSIONAL FLOWS OF AN IMPERFECT DIATOMIC GAS 

By A. J. Eggehs, Jr. 


SUMMARY 

With the assumptions that Berlhelot’s equation of state 
accounts for molecular size and intermolecular force effects, and 
that changes in the vibrational heat capacities are given by a 
Planck term, expressions are developed for analyzing one- 
dimensional flows of a diatomic gas. 

The- special cases of flow through normal and oblique shocks 
in free air at sea level are investigated. It is found that up to a 
Mach number of 10 the pressure ratio across a normal shock 
differs by less than 6 percent from its ideal gas value; whereas 
at Mach numbers above 4 the temperature rise is considerably 
below and hence the density rise is well above that predicted 
assuming ideal gas behavior. It is further shown that only the 
caloric imperfection in air has an appreciable effect on the 
pressures developed in the shock process considered. The effects 
of gaseous imperfections on oblique shock flows are studied from 
the standpoint of their influence on the lift and pressure drag 
of a flat plate operating at Mach numbers of 10 and %0. The 
influence is found to be small. 

INTRODUCTION 

A -wide variety of problems in compressible flow has been 
solved on the assumption that air behaves as an ideal dia- 
tomic gas. This assumption is justified, provided the 
pressure and temperature range of interest is small and near 
atmospheric. It is an experimental fact, however, that 
when air is subjected to large changes in state at pressures 
and temperatures far removed from atmospheric, it ceases 
to obey the simple gas law and exhibits other properties not 
characteristic of an ideal gas. Consequently, flow processes 
in which air is subjected to these extreme conditions can be 
expected to depart from perfect gas behavior. It is known 
that such flows will be encountered in hypersonic wind 
tunnels and by aircraft flying at high supersonic airspeeds; 
hence, the nature and extent of this departure have become 
important considerations in aerodynamics. 

Classical theories and experiments have shown that three 
properties of a real gas first cause it to exhibit characteristics 
unlike those of an ideal gas. These properties may be 
classified as thermal and caloric imperfections. Thermal 
imperfections in the form of intermolecular forces and 
molecular size effects are significantly manifest at low 
temperatures and high pressures. Changes in the vibrational 
heat capacities become an important caloric imperfection at 
relatively high temperatures. Circumstances under which 
effects of molecular dissociation and/or electronic excitation 
become important- (e. g., temperatures appreciably above 


5000° R) may be neglected for the present. Insofar as gases 
in equilibrium are concerned, it is usually sufficient to ac- 
count for intermolecular force and molecular size effects with 
additional terms in the equation of state. Similarly, changes 
in the vibrational heat capacities of the molecules may be 
accounted for with a function of temperature in the expres- 
sions for the specific heats. 

Tsien (reference 1) investigated the effects of gaseous 
imperfections on air flows using Tan der Waals’ state equa- 
tion. Approximate solutions to the one-dimensional isen- 
tropic and normal shock equations were obtained. (Tsien 
points out two very early papers of limited scope on the 
subject by A. Busemann and W. J. Walker.) The Joule- 
Thomson effect was neglected in Tsien’s analysis, 1 however, 
thus introducing some error. Donaldson (reference 2) 
observed this error and found that the differential equations 
of motion could be integrated to yield one-dimensional 
isentropic flow equations, exact to the accuracy of Tan der 
Waals’ equation. A comparison of results obtained with 
these equations and the flow equations for a perfect gas 
indicated that appreciable error would be made if the latter 
expressions were used to predict aerodynamic phenomena 
involving- high temperatures or high pressures. 

A matter considered by Eantrowitz (reference 3) is the 
inability of a gas composed of polyatomic molecules to in- 
stantaneously adjust its internal energy to temperature 
changes at high temperatures. This time lag in equilibrium 
partition of energy occurs primarily in the vibrational 
energies of the molecules which, unlike the translational and 
rotational energies, require many collisions per molecule to 
become fully adjusted to a new temperature. As pointed 
out by Eantrowitz, for rapidly accelerating flows the result- 
ing heat-capacity lag causes an entropy increase and hence a 
total-head decrease in a gas. The net effect may be made 
negligible, however, by decreasing the rate of change of 
temperature of the fluid with time. In the special case of 
flow through a normal shock wave, Bethe and Teller (refer- 
ence 4) found that heat-capacity-lag effects caused the thick- 
ness of the shock to increase. (The wave is considered to 
extend over that portion of the fluid in which equilibrium 
partition of energy does not prevail.) As the Mach number 
upstream of the wave increases, however, the shock tends 
rapidly toward a sharp discontinuity. It may also be 
observed that, although heat-capacity lag influences flow 
within the shock, it- does not alter the magnitude of the 

1 This error is reported to have been, subsequently corrected; however, the author is una- 
ware of the publication. 
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entropy increase across the shock, as this is independent of 
the mechanism by which available energy is dissipated within 
the wave. These observations are important for they lend 
support to the simplifying assumptions that moderately 
accelerating flows may be treated as isentropic except in the 
presence of a shock wave, in which case the wave may be 
treated as a discontinuity and the properties on both sides 
calculated as if there were no heat-capacity lag. 

In the present paper, one-dimensional flow equations are 
obtained with the aid of Berthelot’s equation of state. This 
equation, rather than Van der Waals’, is employed in order 
that somewhat better estimates of intermolecular force 
effects may be obtained. It is assumed that caloric imperfec- 
tions may be accounted for with a Planck term in the ex- 
pressions for the specific heats. The analysis covers both 
isentropic and plane-shock flows. Within the limitations . of 
the assumptions, both exact and approximate flow equations 
are developed. The approximate expressions provide explicit 
solutions for all the important flow parameters, and are 
designed to utilize data available on the theoretical behavior 
of an ideal diatomic gas. (See, e. g., reference 5.) 

SYMBOLS 

a local speed of sound, feet per second 
A nozzle cross-sectional area, square feet 

b molecular size constant, cubic feet per slug 

c intermolecular force constant, °R, feet to the fifth per 

slug, second squared 

c„ specific heat at constant pressure, foot-pounds per 
slug, °R 

e, specific heat at constant volume, foot-pounds per 
slug, °R 

e base of natural logarithms, 2.718. . . 

In logarithm to base e 

M Mach number (ratio of local velocity to local velocity 
of sound) 

P absolute pressure, pounds per square foot 
R gas constant ^1715 for foot-pounds per 

slug °R 

T absolute temperature, °R 

u internal energy, foot-pounds per slug 

V local velocity, feet per second 

v specific volume cubic feet per slug 

w external work performed, foot-pounds per slug 

Y ratio of specific heats ( c p jc v ) 

p mass density, slugs per cubic foot 

6 molecular vibrational energy constant, °R 


ANALYSIS 


It will be assumed throughout the analysis that the gas 
exhibits no heat-capacity lag. Shock-free flows are con- 
sidered isentropic, and flow through plane shock waves is 
assumed to be adiabatic. Justification for the first two 
assumptions is discussed in the Introduction. The extent 
to which plane shock processes, particularly at high Mach 
numbers, deviate from adiabatic behavior is left for future 
investigation. 

Berthelot’s equation of state is employed in two forms, 
depending upon the desired range of applicability. In the 
development of equations for investigating flows over a wide 
range of Mach numbers, temperatures, and pressures, the 
state equation is used in its exact form 


pRT cp 
1-bp T 


where b is the molecular size constant and c is the intermolec- 
ular force constant. It will be noted that equation (1) 
differs from Van der Waals’ equation by a factor of 1/Tin the 
intermolecular force term. The introduction of this factor 
yields a variation of intermolecular force with temperature 
which is in accordance with experiment. (See, e. g., refer- 
ence 6.) Thus, as pointed out in reference 7, close agree- 
ment between theoretical and experimental values of com- 
pressibility of gases is obtained. 

In the derivation of approximate flow equations an expres- 
sion of the form 

P~ P RT(l+bp--$p) ( 2 ) 

is employed. This equation is, of course, a first-order approx- 
imation to equation (1). 


FLOW OF A DIATOMIC GAS OBEYING BEETHELOT’S EQUATION OF STATE 

Isentropic flow of a Berthelot gas may be studied with the 
aids of the differential isentropic expansion equation, 


where 



and the Joule-Thomson effect is given by 



( 4 ) 


SUBSCRIPTS 

0 stagnation conditions 

1 conditions upstream of shock wave 

2 conditions downstream of shock wave 

i ideal gas quantities 

SUPERSCRIPTS 


Combining these expressions and substituting the value of P 
from equation (1) yields 

c„dT ijsdp dp— 0 (5) 

Now the differential expression for c, is 


quantities at the nozzle throat 
quantum-mechanical functions 



P 
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which, upon substituting from equation (1), may be inte- 
grated to give 

( 6 ) 


c —c 

W — w v I /Jt 2 


where c/ is a function which describes the quantum-mechan- ' 
icai variations of c, with temperature. The second term on 
the right of equation (6) represents the effects of gaseous 
imperfections on c„. 

The function chosen for c,' is determined by the molecular 
structure of the gas under consideration and the temperature 
range over which accurate predictions of c/ are desired. Tor 
aerodynamic purposes, diatomic gases are of primary interest. 
The important temperature range extends from liquefaction 
temperatures to several thousand degrees Eankine. A rela- 
tively simple function may be written for c/ in this case, as 
the number of translational and rotational degrees of freedom 
is constant, and only the variation with temperature of the 
vibrational heat capacity need be considered. This function 
is 

( / flV e w/T> ) 

c/=Cj t jl + (7 i 1) ^_ e c9/r)j2j (7) 

The second term in the brackets, essentially a Planck term, 
accounts for the vibrational contribution to the specific 
heat at constant volume. 2 The assumption is that the mole- 
cules of the gas behave like linear harmonic oscillators insofar 
as the vibrational degrees of freedom are concerned. (See 
reference 8.) 

An expression governing isentropic expansion of an imper- 
fect diatomic gas may now be obtained by substituting 
equations (6) and (7) into equation (5) and integrating from 
stagnation to static states. Since 

Tdp — 2 pdT 


: ( T *)“■ 


T* 


there results the relation 

, Tn ( T \± Pin M1-&P) - ( P _ A> \ , Wn 1 

c H ln (j-J+Pdn p (i_ 6po) T z)+Rln 


R 


fe e« tT > T , 
[_T 


( 8 ) 


In order to deter min e the Mach number of a stream, it is 
necessary to find the velocity of flow and speed of sound in 
the stream. These quantities may be found by employing 
the one- dim ensional energy equation, 

9» + «p»)+iw-@)y(i)+(^)ir+9(f)+w=o 

(9) 

Substituting equations (1), (4), (6), and (7) into equation (9) 
and integrating from stagnation to static temperature and 
density yields for the velocity 


V>=2{c h (T 0 -T)+^ 


Ro 


e com] 


+2c 


r 0 




(10) 


The corresponding speed of sound is determined by substitut- 
ing equations (1), (5), (6), and (7) into the general equation 

2 _dP_fdP\ ,/aP\ dT 

a -dp-\bp) T + \j!>T), dp 
The resulting expression is 


RT 2 c P , 

p 2 T 

1± 

R 1 
p(l-&p)_ 

2 

a (1 — M 2 

t— 1 
1 

+ 

iH 


* , 2cp] 

' [l_ e («mj2Tj5y2i 

n 

u 


2 This is a common method of accounting for the variation with temperature of the vibra- 
tional heat capacities. It has been adopted by Donaldson and others for imperfect gas 
studies. 

Combining equations (10) and (11) yields the following equation for the Mach number 

(v— 1 ) (jj Tp + TZT*) , l 2 cp 


(U) 


M i =2< 




(1 — 6 p) 2 RT- 


( 12 ) 


The specific beats and the ratio of specific heats are readily obtainable for a diatomic gas obeying Berthelot s state equation. 
The specific heat at constant volume c„ is found by substituting equation (7) into equation (6), thus yielding 


r. , , « w/r) , 2 cp r 

c,—c u l + (yt— l)Ujy [i — gW/^pT" RT 2 ) 


(13) 


The specific heat at constant pressure c„ is obtained by substituting equations (1) and (13) into the reciprocity relation 

and setting 


(dP/azy 

* -Cs 1 (bPfin ;), 


c Vi — c Pi R 


(14) 
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The resulting expression for c p is 


Cp — 


1+ 



- 

/2 —bp ,1 cp \~ 

*> 


/0Y e (sm , 2cp 

, , \1 -bp 2 RT 2 ) 

L 

\ r i ) 

\TJ 1 RT 2 

" r 1 2cp 



{ 

(1— 6 P ) 2 RT 2 



( 15 ) 


The ratio of specific heats y follows directly, of course, from equations (13) and (15) and may be written 


y=yi 



f 

(2— bp , 1 cp \ 

■v 

j-A* 1 . 

/0Y e w/r) , 2 cp 

\l-bp 1 2 RT 2 ) 


1 ' 
7 i 

\TJ [1— e w/T) ] 2 1 RT 2 

1 1 2 c P 



L 

(1 — 6p) 2 RT 2 



,,, * mT) ,2 cp\ 

L-K7i l) \\T)[l-e^ T) f + RT 2 ] 


(10) 


A matter of special interest in the study of channel flows 
is the variation of Mach number with cross-sectiona] area of 
the channel. The relation between these two quantities is 
most conveniently determined from the ratio of cross- 
sectional area at the sonic velocity station (the station of 
minimum area for any gas, as pointed out by Tsien) to 
cross-sectional area at a station of arbitrary Mach number. 
The continuity equation yields this ratio in the form 




(17) 


With the aid of this equation and equations (1), (8), (11), 
and (12), the area ratio corresponding to a particular Mach 
number may be determined from a knowledge of the ideal 
gas specific heats, characteristic constants of the gas (i. e., 
b, c, and 0), and, for example, the stagnation conditions. In 
general, however, this computation cannot be earned out 
entirely analytically. For example, it is seen from equation 

(8) that none of the variables can be explicitly determined; 
hence, a final solution for any single variable must be ob- 
tained graphically or by equivalent means. 

In order to relate the fluid properties on the two sides of 
a normal shock wave, it is convenient to integrate equation 

(9) , the energy equation, across the wave. The resulting 
expression is 


(§-§> 


i 


(18) 


Continuity of flow and conservation of momentum must also 
be satisfied across the wave. These requirements may be 
expressed analytically in .the familiar forms, respectively, 


and 


piVi — P2V2 


P 1 V 1 2 -p 2 F 2 2 =P 2 -P, 


(19) 

( 20 ) 


Equations (18) through (20) provide the additional informa- 
tion necessary to determine the flow through a normal shock 
wave. Here again, however, final solutions for certain of 


the unknown flow parameters are most conveniently obtained 
graphically. 

APPROXIMATE FLOWS OF AN IMPERFECT DIATOMIC GAS 

This part of the study is concerned with first-order devia- 
tions of one-dimensional gas flows from the behavior of an 
ideal gas. With this restriction, explicit solutions can be 
obtained for all the important flow parameters. Equation 
(2) will, of course, be employed as the equation of state. A 
simplification is also allowed in the temperature function for 
the vibrational contribution to the specific heats, since the 
specific heats are restricted to vary only in the first order' 
from their ideal gas values. In this case only large values 
of 9/T are considered and the Planck term may be approxi- 
mated by 


(8V_e^_ / 0 V wn 

\T/ [l-e WT >Y-\Tj e 


Equation (8) may now be reduced to the form 
1 


-l]+ 


^[(to) ’ (:)W: o + mv [(T/To) 2 

(ijr+l) c~ w3 V-(j,+ 1 ) (21) 

This equation can be readily solved for p/p 0 in terms of T/T a) 
retaining only terms of the proper order, thus yielding 

I 1 3—2 Vi 

H-KID* ~ i ]~#&[(£)’ _i ]~ 


Po 


(t^+i)«-«™+[(3t) (p)+l]e- wr «>™ r) } (22) 

By solving for TjT 0 in a similar manner, there is -.obtained 

To = (£) ' [ 1+b_1) { 6p °[(^) _ 0 + ^’[(po) ' -1 ] + 


( 23 ) 
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A relation between the pressure, and density ratios for isen- 
tropic expansion is obtained by comb inin g equation (2) with 
equation (23). The resulting expression is 


erCD f 1 0 : 


-(®T U > 


[ttxr+i] 


.-(£®r 


}] 


(24) 


Equations (22) through (24) are useful for determining the 
nature and extent of small departures from ideal gas behavior 
of a diatomic gas undergoing isentropic changes in state. 

Isentropic relations will now be obtained between pres- 
sure. temperature and density ratios, and the stream Mach 
number. In order to do this, it is convenient to reduce 
equation (12) to a form consistent with the first-order 
approximations. Performing this operation yields 




'« ( P/Pa\ - c Pa ( -i p/po ^ , 

'ToV p/po / btA t;tJ + 


g-tf/Vo) C T 0 m 


PoBZ 


JL ^g-W/ro) g-C«/r IJ )(i’Q/r)j ^ 

To 


where 


) 


r i =(7r-3 Tt +3) 


(25) 


An expression for the temperature ratio in terms of Mach 
number, stagnation temperature and density, and gas con- 
stants is obtained by substituting equations (2), (22), and 
(24) into equation (2a) and solving; hence 

Y 0 = (%)i [l+5i (AO&po+CUJi) j$k+ 

where 


and 


- 2 s (£>? D-GDn-* 7 -w L 1 - ©n 


(27) 


where (_T/T 0 )t, (p/p 0 ) { , and (P/P 0 ) t are unique functions of the Mach number M. (These functions and all other ideal gas 
expressions employed in the analysis are included in Appendix A.) The density and pressure ratios are obtained by sub- 
stituting equations (23) and (24), respectively into equation (26), thus yielding 


where 


and 


Tf— 1 


RTj 

where 


P 3 (M> 


B l (M) 


S,(M) = l-(^) 

\Po/j Ti — 1 

d 2 (Jt,m) (i +^°)+ 


(28) 


(29) 


Similarly 


rO,[ 1+ftMt '* +Gm w 

D >(r;- W )(r t )'-'™] 


( 30 ) 


and 


7i— 1 

( 1 + t)+ 


[(t) + T°] e ~ w/ro> 1 WT, *“ 11 + 


Y (31) 


7 £ 

It may be noted that equations (26), (28), and (30) can be 
rearranged to give- the Mach number in terms of TjT 0l p/po, 
and P/Po, respectively. 

Expressions for the specific heats and the ratio of specific 
heats are derived by reducing equations (13), (15), and (16) 
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to conform to the first-order deviation theory. Performing 


this operation yields 



e.«e„{l+(Y«-l) 

2cpl\ 

RT* J) 

(32) 


6c P "|) 
~RT* J) 

(33) 

and 



V=Vi{ 1 -(^)[(7i- 1) + frr 

-3)|^]| 

(34) 


The local speed of sound a can be obtained directly from 


It is observed that c P) c Vi and y are independent of molecular 
size within the limitations of the theory. This fact appre- 
ciably simplifies the determination of the characteristic gas 
constants c and b. 

The calculation of cross-sectional area of a nozzle for 
isentropic expansion of an approximate Berthelot gas can 
now be made. For this purpose it is convenient to write 
equation (17) as follows: 

S-COGDGO* 

ion (11) by retaining only first-order terms, thus 

(35) 


Hence, substituting equations (26), (28), and (35) into equation (17) in the above form yields for A* /A, 



l+B.mbpo+G^M) 


Cpo 


BTa‘ 




(36) 


where 




y (37) 


-M/T 0 ) 


{(#),+$ (£)'($0 


and {A* l A) { is the explicit function of Mach number given in Appendix A. 

Flow through a normal shock wave will now be considered. The temperature, density, and pressure ratios across a 
normal shock wave may be related to the ideal gas expressions for these quantities by simplifying and rewriting equation 
(18) as follows: 


(i-f )+Tf‘ [(^)~i]+ ( i -Ps)+» 


(38) 


Substituting equations (2), (19), and (20) into equation (38), dividing through by equation (35), and solving for the tempera- 
ture ratio yields the expression: ; 


h 

2V 


<t)i 1+ ^T^ — [B.(M,)i Pl +C,(M l) A(£.M,)(£) 

\ Vi — l 




(39) 


where 




and 


-(£> {(ID.-©. ,+^i [>-(§), ] +i j 

{3(g) [$&-i]+£i [i-(fo ,])■ -©hi:-© ij-i}-©,©; 

(^’ ■“-)= f C^)<{ i -C-r«-u(£) [i -(^) j]+ ^-0 (^) (g), [i -(g) J 


( 40 ) 
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In order to obtain the relation for the density ratio p 2 /pi, 
it is convenient to write equation (39) in the form 

t=d), ' 1+ -> (41) 


where e T includes all the first-order correction terms in 
equation (40). Combining this expression and equations 
(2), (19), and (35) with equation (38) yields the relation 


H£), l 1 -®', W [W0*»+<WK> sfr 


where 


Ti+1 


*4!^ KS)r i ] +s [(twit 1 ]} 


and 


D e (4r> MiV- fl + 


(42) 


(43) 


As with the temperature ratio, it will be convenient to write 
equation (42) in the form 


?“(?)< Cl+ *) (44) 

Pi \Pi/» 

where e p is the sum of all the first-order correction terms. 
The pressure ratio across a normal shock wave is now ob- 
tained by combining equations (19), (20), (35), and (44) 
to provide the relation 


Stagnation conditions downstream of the wave may be 
related to stagnation conditions upstream of the wave with 
the three identities 

P o 2 /PA / Pl/P °A 
Po~\Pi) \P 2 IPj 

(49) 

a >2 / pA f pi/poA 
Po t VPi/ VP2/P0J 

(50) 

and 

To 2 /TA (TJTA 
T a ~\Ti) \T 2 /tJ 

(51) 


in which all terms of the right-hand members are given in 
previously developed expressions. 

It is important to note that the normal shock expressions 
are applicable to the study of oblique shock flows. For 
example, equations (39) through (46) may be employed 
directly by substituting M\ sin jS for Mi wherever it occurs. 
In order to dete rmin e Mi, however, it is convenient to have 
a relation between the wave angle j3 and the stream deflection 
angle 5. This relation may be obtained by combining the 
momentum and continuity equations for flow through an 
oblique shock (reference 9) to yield the expression 

P» , _pjYI ( i \ 

Tt Pi Vcot P cot 5+1/ 


which may be written 

cot 5=tan/3 jjy*ZWi a (l — i?) (Pg/pj — 1 ~ 1 
Mi can then be calculated from the expression 

fpi) 

Mi \Tj t \piJi sin 03—5) 


(52) 


(53) 


in 



fe+ 17 )— (45) 
(£) (46) 


The additional information needed to define the flow 
through a normal shock wave is the Mach number Mi on 
the downstream side of the wave. This Mach number may 
be obtained with the continuity equation 

■Ml = ( flA ( £A 

Mi \Pi) \&2/ 


by substituting in equations (35), (41), and (44), thus yield- 
ing 



(47) 


where 


_ «r,i fi /pA 1 1 A c Pi 1~] (Pa/Pi) t ~| 

t “—~' 2~ rbpl L 1 - W< _Tt* BT? L OTTOAJ 


+ 


hi-iY f & 


27 * 


(xJ { (Mil e ~ w/T ° aTi/T2) ~ 1! _ 1 1 (48) 


which is obtained from the continuity equation. In this 
equation the ideal gas values for the temperature aDd density 
ratios are those corresponding to the initial Mach number 
Mi and the wave angle j3. "When M 2 from equation (53) is 
known, the remaining normal shock expressions, equations 
(49) through (51), can also be applied to oblique shock flows. 

It is observed that, in general, all expressions obtained in 
the preceding approximate analysis differ from the corre- 
sponding ideal gas equations only by small correction factors 
containing the Mach number, temperature and density, and 
constants of the gas. As pointed out in the Introduction, 
this simp lifi es the investigation of one-dimensional imperfect 
gas flows, provided effects of the imperfections are not large. 

DISCUSSION AND CONCLUSIONS 

The complex manner in which gaseous imperfections influ- 
ence general one-dimensional-flow processes precludes the 
detailed consideration, here, of more than two flows of aero- 
dynamic interest. Before discussing these flows, it may be 
worth while to establish more definitely the conditions under 
which imperfections in air can alter its dynamical behavior 
from that predicted on the assumption that ideal gas laws 
are obeyed. These conditions are easily deduced from the 
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Figuee 1.— Variation of pressure ratio across a normal shock wave with Mach number, 
(Initial state— N AC A standard atmosphere at sea level.) 


equations employed in the analysis and a knowledge of the 
characteristic gas constants . 3 They may be summarized as 
follows: 

1. The occurrence of temperatures in excess of 800° R, in 
which case the specific heats change appreciably. 

2 . The occurrence of densities in excess of 10 times sea 
level atmospheric, in which case the molecules occupy an 
appreciable fraction of the available volume. 

3. The occurrence of sufficiently high pressures and low 
temperatures to cause the intermolecular forces to be 
appreciable. 

One additional observation should be made concerning the 
integrated effect of ..these phenomena on a particular flow 
process. This is that the several imperfections frequently 
counteract each other (note, e. g., the state equation and 
equations (24), (26), and (27)), and thus alter certain flow 
parameters only slightly from their ideal gas values. 

The special cases of flow through normal and oblique 
shocks in free air (NACA standard atmospheric conditions 
at sea level) will now be considered. The pressure and 
temperature ratios across a normal shock have been calcu- 
lated for Mach numbers from 1 to 10 using equations ( 1 ), 
( 11 ), (18), (19), and (20). In this case the terms containing 
b and c are extremely small (it is easily shown that they could, 
in fact, be neglected), and it is sufficiently accurate to use 
the values for these constants given in Appendix B. 

The pressure rise so determined is shown in figure 1 as 
a function of Mach number, and it is seen that it differs by 
less than 6 percent from the ideal gas value 4 up to the maxi- 

3 The values chosen for the characteristic gas constants, and the method by which they 
were determined, are included in Appendix B. 

4 All ideal gas data are obtained from the tables of reference 5. 


mum Mach number of 10 . This behavior might be antici- 
pated from the ideal gas equation for the pressure ratio 
(Appendix A) which is relatively insensitive to variations 
in 7 at high Mach numbers. The corresponding temperature 
rise across a normal shock is illustrated in figure 2 . Here it 
is observed that the effect of increase in specific heats appre- 
ciably reduces the temperature ratio below the ideal gas 
values at the higher Mach numbers . 5 It may be deduced 
from figures 1 and 2 , that the density rise will he considerably 
above the ideal gas value at the higher Maeli numbers. 
The accuracy of these results will decrease at Mach numbers 
in excess of 7, for under these circumstances the temperatures 
downstream of the wave appreciably 'exceed 5000° R which 
is sufficiently high to cause molecular dissociation and thus 
alter the flow. It is indicated in reference 4 , 5 however, 
that up to a Mach number of 10 , this decrease in accuracy is 
significantly manifest only in the temperature and density 
ratios. (For example, at M= 10 , the effect of dissociation is 
to further decrease the temperature ratio from 17 to about 
14.5.) 

In view.. of the fact- that the changes with temperature of 
the specific heats of air do not strongly influence the pressure 
rise across a normal shock, up to Mach numbers of 10, it 
follows that corresponding pressure effects in flow through 
oblique shocks wilL also be small, provided the Mach num- 
ber of the component of velocity normal to the wave is below 
10 . This suggests, within the limitations observed, that 
the effect of variable specific heats on those aerodynamic 



Floats 2.— Variation of temperature ratio across a normal shock wave with Mach number. 

(Initial state— NACA standard atmosphere at sea level.) 

s This result is consistent with Wood’s findings concerning surface temperatures al high 
Mach numbers. (See reference 10.) 

8 In this reference, calculations similar to those described above have been carried out 
including the effects of molecular dissociation but neglecting the oilocts of thermal imperfec- 
tions, The results of this reference and the present report arc found, where comparable, to 
be in good agreement. 
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characteristics essentially independent of viscosity may be 
small for some shapes. (Epstein’s findings in reference 11. 
concerning hypersonic flows about bodies, pertain to much 
higher Mach numbers than those considered here.) The 
extent to which this is correct has been investigated for the 
simplest shape, a flat plate. Pressures on the upper surface 
of the plate were determined, assumi n g Prandtl-Meyer 
flow about the leading edge, 7 and pressures on the lower 
surface were found with equations (45) (replacing Mi with 
Mi sin jS) and (52) . s These data and the corresponding lift 
and pressure drag coefficients of the plate were calculated 
for Mach numbers of 10 and 20, and for angles of attack 
from 0° to 24° and 13°, respectively. A comparison of the 
coefficients with those obtained assuming ideal gas flow 

7 The assumption of Prandtl-Meyer flow- trill yield increasingly inaccurate upper-surface 
pressure as t he local stream temperature approaches absolute zero. In this case the pressures 
are so low, however, that they make a negligible contribution to the forces on the plate. 

5 This method of calculating lower-surface pressures requires that the chord of the plate be 
large in comparison to the region of the leading edge contained within the shock wave. 


over both the lower and upper surfaces showed that the 
effect of variable specific heats is very small. (For example, 
the force coefficients are smaller in the case of the imperfect 
gas but never by more than 3 percent, even at the highest 
angles of attack.) This result is not justification for as- 
suming ideal gas behavior to calculate the high Mach 
number aerodynamic characteristics of arbitrary shapes, 
other effects such as specific heat lag being neglected; but 
it does indicate that the assumption is justifiable for certain 
shapes on which the important pressure forces are obtained 
through oblique or normal shock compression. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffet Field, Calif., Dec. 30, 1948. 
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APPENDIX A 
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IDEAL GAS EQUATIONS 

The following are ideal gas equations used in the analysis 
of approximate flows of an imperfect diatomic gas: 

ISENTBOPIC FLOW EQUATIONS 
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APPENDIX B 


DETERMINATION OF CHARACTERISTIC GAS CONSTANTS 

FOR AIR 

Numerical values for the characteristic gas constants b 
and c may be determined from a knowledge of the critical 
pressure and temperature of a gas. Values determined in 
this manner are, however, generally less accurate in a par- 
ticular range of temperatures and pressures than is desired. 
Consequently, it is often advantageous to fix b and c such 
that theoretical and experimental values of the primary 
parameters c pi y and compressibility (i. e., P/p) are in good 
agreement in the range of interest. This method is the 
more suitable of the two for most aerodynamic studies as 
the variations in pressure and temperature are limited; 
hence, it has been employed in this paper to determine b 
and c, as well as 0, for air. It is sufficient for the illustrative 
applications presented here to determine the former con- 
stants only for the approximate flow equations. Separate 
values of 0 are found for the exact and approximate Planck 
terms. 

As pointed out previously (note the development of approx- 
imate expressions for c„, c v , and 7), it is a simple matter to 
obtain the characteristic gas constants appearing in the ap- 
proximate flow equations. In this example, c and 0 were 
chosen first to give agreement between experimental and 
theoretical values of c P and 7, the latter values being, calcu- 
lated from equations (33) and (34), respectively. Data on 
the variation of these quantities with pressure and tempera- 
ture at high temperatures were obtained from references 12 
through 19.® A comparison of these data with theory, setting 

8 These data have been correlated by the Research Department of the Pitts burgh-Des 
Moines Steel Company, and the results of this correlation are presented here. 
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c— 2.25X10 8 °R, ft 5 /slug, sec 2 and 0=5800° It, is shown in 
figures 3 and 4 for pressures of 0 and 144 X10 9 pounds per 
square foot absolute. The agreement is observed to be good 
up to temperatures of 3000° R. Using equations (15) and 
(16) to calculate c„ and 7, respectively, and a value of 
0=5500° R, it is seen that at zero pressure excellent agree- 
ment with the correlated data is obtained up to 5000° R. 
The approximate theory was then checked at low' tempera- 
tures with experimental values of 7 given in reference 20. 
At temperatures above liquefaction and pressures up to 25 
atmospheres, the difference between theoretical and experi- 
mental values of 7 was found not to exceed 3 percent. 
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Figure 3.— Variation of C, with temperature at absolute pressures of 0 and 144X10* pounds 
per square foot absolute. 
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Figure 4.— Variation of r with temperature at absolute pressures of 0 and 144X10 3 , pounds per 
square foot absolute. 


The molecular size constant 6 was chosen to yield good 
agreement between theoretical and experimental values of 
compressibility of air. The theoretical values were calcu- 
lated from equation (2), and the experimental data were 
obtained from reference 21. In general, it was found that 
using 6=0.485 ft 3 /slug and the previously determined value 
of e, equation (2) predicts the compressibility of air to within 
2-percent error at pressures from 0 to 150 atmospheres, and 
temperatures from liquefaction to 850° R. (The deviation 
from ideal gas behavior is small above this temperature at 
the pressures indicated.) 
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